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Abstract. We calculate the bispectrum of the Cosmic Microwave Background 
(CMB) temperature anisotropics induced by the second-order fluctuations in the 
Boltzmann equation. In this paper, which is one of a series of papers on the numerical 
calculation of the bispectrum from the second-order fluctuations, we consider the terms 
that are products of the first-order perturbations, and leave intrinsically second-order 
terms and perturbations in the recombination history to the subsequent papers. We 
show that the bispectrum has the maximum signal in the squeezed triangles, similar 
to the local-type primordial bispectrum, as both types generate non-linearities via 
products of the first-order terms in position space. However, detailed calculations show 
that their shapes are sufficiently different: the cross-correlation coefficient reaches 0.5 
at the maximum multipole of Imax 200, and then weakens to 0.3 at Imax ~ 2000. 
The differences in shape arise from (i) the way the acoustic oscillations affect the 
bispectrum, and (ii) the second-order effects not being scale- invariant. This implies 
that the contamination of the primordial bispectrum due to the second-order effects 
(from the products of the first-order terms) is small. The expected signal-to-noise 
ratio of the products of the first-order terms is ~ 0.4 at Imax ^ 2000 for a full-sky, 
cosmic variance limited experiment. We therefore conclude that the products of the 
first-order terms may be safely ignored in the analysis of the future CMB experiments. 
The expected contamination of the local-form /jvl is /^£°' ^ 0.9 at Imax ~ 200, and 

/^Y' - 0.5 at Irnax ^ 2000. 
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1. Introduction 

Primordial non-Gaussianity is now recognized as a powerful probe of the details of the 
physics of inflation [1], as detection of large primordial non-Gaussianity would rule out 
all classes of inflation models that satisfy the following four conditions simultaneously: 
single-field, canonical kinetic term, slow-roll, and initially vacuum state. 

However, the extraction of the primordial non-Gaussianity may not be so simple, 
as there are various non-primordial effects that can also generate non-Gaussianity. Any 
non-linearities can make initially Gaussian perturbations non-Gaussian. 

The angular bispectrum, Bi-^i^i^, the harmonic transform of the angular three-point 
function, of the Cosmic Microwave Background (CMB) is often used to measure non- 
Gaussianity (see, e.g., [2], for a review). Departures from any of the above conditions 
(single-field, canonical kinetic term, slow-roll, and initial vacuum state) would result in 
detectable non-Gaussian signals in specific triangle configurations of the bispectrum. 

When we consider the effects of various non-primordial sources of non-Gaussianity 
on the extraction of the primordial signals, we must specify of which primordial non- 
Gaussianity we study the contamination from the non-primordial sources. Multiple- 
field models, non-canonical kinetic terms, and initially excited states can produce large 
signals in the squeezed triangles {li <^ h ^ h) [2], the equilateral triangles (/i = I2 = h) 
[1], and the flattened/folded triangles (/i ^ I2 ^ hf^) [SIE], respectively. 

Throughout this paper we shall study the contamination of the squeezed triangles, 
parametrized in the form of the so-called local form of the bispectrum, which results from 
the primordial curvature perturbation (in comoving gauge) in position space, C(x)) given 
by C(x) = (L{'x) + ^fNL(1{x.), where Cl is a Gaussian perturbation, and /atl characterizes 
the amplitude of the local-type non-Gaussianity. Our sign convention is such that the 
temperature anisotropy in the Sachs-Wolfe limit at the first-order in perturbations is 
given by AT'-^-'/T = (1/5)^^^-*. The simplest class of inflation models satisfying all of the 
four conditions (single-field, canonical kinetic term, slow-roll, and initial vacuum state) 
produce very small non-Gaussian signals: f^L ~ 10~^ at the horizon crossing during 
inflation [71 [8], whereas the best limit from the WMAP 5-year data with the optimal 
bispectrum estimator is /at^ = 38 ± 21 (68% CL) [9]. How much would non- primordial 
contributions account for the measured value of /atl? 

The CMB bispectrum from the local-type primordial non-Gaussianity with the 
linear radiative transfer has been given in [10], and that arising from non-linearity 
in gravity has been considered in [11]; however, non-linearities exist also in the 
perturbations in the photon-baryon fluid, i.e., non-linearities in the Boltzmann equation 

In this paper we calculate the CMB bispectrum, taking into account the second- 
order perturbations in the Boltzmann equation. We shall include the second-order terms 
that are products of the first-order perturbations, and ignore the intrinsically second- 
order terms (some of them have been considered in [T5l [T6l [T7]). or the effects of the 
perturbed recombination [T8| [T9 | [20] . The calculations that also include the intrinsically 
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second-order terms and the perturbed recombination will be presented elsewhere. 



2. CMB Bispectrum From Second-order Perturbations 

2.1. Definitions 

We expand the temperature fluctuation into the linear (first-order) part and the second- 
order part as 

AT(n) _ AT(^(n) AT(^)(n) 

rp T ^ T ' ' ' ' 

The spherical harmonic coefficients of temperature anisotropy, ai^ = J d'^iiYi^{h)AT{h) 
are therefore expanded as 

a,m = + + • • • • (2) 

(2) 

How do we calculate the second-order part, a^^? This can be calculated by expanding 
the Boltzmann equation up to the second order in perturbations |12j . 

To expand the Boltzmann equation up to the second order in perturbations, we 
first expand the distribution function, 

n -1 



/(x,p,n,?7) = 2 



(3) 



up to the second order in perturbations: = O*-^^ + 0^^''/2 + . . . , and accordingly 
/ = /(0) + /(i) + ;(2)/2 + .... 

We compute the fractional perturbation in photon's energy density at the i-th order 
in perturbations, A^*^ by multiplying /*^*) by p, and integrating over p'^dp: 

I dpp^fi^) ■ ^ ^ 

At the linear order, we recover the usual relation between the linear fractional 
temperature fiuctuation, B*^^-* = AT^-^^T, and the linear fractional energy density 
perturbation, A^^) = 6p^^^/p^, i.e., A^^) = 4e(^). 
At the second order we have 

A(2) = 4e(2) + i6[e«]2, (5) 

which is related to the second-order temperature fiuctuation as 

= 1 (A(^) - (A(^))) - 1 ([e«]^ - ([e«]^)) 

= l(0(2)_(0(2)) + [0(l)]2_^[@(l)]2))^ (6) 

where we have subtracted the average of the temperature fiuctuation so that the average 
of AT^'^^T vanishes. 
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We compute from AT^^^^T using 



2 

I'm' l"m" 



where we define 



a^2 ^IJ d'hY,Uh){A^'\h) - (A(2)(n))), (8) 



(2/i + l)(2/2 + l)(2/3 + l) h k h\( h h h 
An I / I mi m2 



(9) 



Here tlie matrix is tlie Wigner 3j' symbol. 

The CMB angular-averaged bispectrum, Bi^i^i^, is related to the ensemble average 

of '^/imi '^/2"l2 '^'3"l3 

Bhhh = ( mg ) ^"'i^i^'a^a^'amg)- (10) 

all m \ / 

This definition guarantees rotational invariance for the bispectrum, and the Wigner 3j 
symbol ensures that the bispectrum must satisfy triangle conditions: — < h < k+lj 
for all permutations of indices, and selection rules: nii + m2 + m-^ = 0. 
The ensemble average is given by 

X ((aS:i,«!la!^i^4L^) - («!:i,«!l)(4i^4L^)) + cyc/zc, (11) 
where cyclic means that we have to sum the cyclic permutations of Eq. ffTTj) for indices 
(1,2, 3) -^(3, 1,2) ^(2, 3,1). 

''Im 



As we assume that a|^^'s are Gaussian random variables, the four-point function of 
a[^]'s in Eq. ( ITT]) is given by the sum of products of all possible pairs. Each pair gives 
the angular power spectrum, Cf. 

= {-irCAv5.mm'. (12) 

We obtain 

\"Zimi "Z2m2 \m'^ I'im'l ' "/2m2 / \"/!jm^ I'im'^ ' 

= (-l)'"^+'"^G,a2[5/l/^^-n^l„^^^^2^^,'^~m2< + (1 2)]. (13) 

Substituting the right hand side of equation (fT3l) for the second term of equation 
( ITTj) . and using h + h + h = even, we obtain the angular averaged bispectrum, 

^/i/2«3 = ^/i/2/3 - 3/«i«2/3(C/,C/2 + cycUc), (14) 
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where we have defined the quantities, 



^ (2/i + l)(2/2 + l)(2/3 + l) / k k h 
'^1^2/3 - A/ 4^ 



(15) 



and 



all m 



2.2. Angular averaged bispectrum from second- order perturbations 

The Boltzmann equation governs the evolution of A^^^ {k, fi,r]) and A'^^)(k, n, 77), where 
fi = k ■ n and n is the direction of propagation of photons. Note that for the hnear 
perturbation there is azimuthal symmetry such that A^^^ depends only on the angle 
between k and n; however, for the second-order perturbation there is no such symmetry. 
The Boltzmann equations in Fourier space are given by 

A«' + 2A;/iA(^) - r'A« = ^«(A;, /i, r/), (17) 
A(2)' + zA;/iA(2) - r'A(2) = S^^\k,h,r]), (18) 

where the primes denote derivatives with respect to the conformal time d/dr], S^^^ and 
are the source functions at the first and the second orders, respectively, and r' is the 
differential optical depth which is defined by using the mean electron number density. 
He, the Thomson scattering cross-section, ax, and the scale factor, a, as 

T = —UgaTa. (19) 

We expand the angular dependence of A*^*^ as 



A!£(k,^) =^Y^ / d'hY,l{h)A^HKn,v), (20) 



and that of the source terms as 



Si'liKv) = / d'hY;^{h)S^'>{k,h,v), (21) 

where i = 1,2. 

The source functions relate the observed ai^s to the primordial curvature 
perturbations in comoving gauge, C(k). The relations contain the linear radiation 
transfer function, gi{k), and the second-order radiation transfer function, F[^' (k), and 
are given by 

d^k 



= I 7^^?/(^)>^/™(k)C(k), (22) 



A-K I f d^k f d^k 



~(2) 



J2 ^L™'(k', k", k)lir^,(k)C(kOC(k"). (23) 



X 

I'm' 
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The linear transfer function is given by 

rvo r ^ Z'? r/^ 1 \ " 

Sg'{k,,i) + S'^i\k,,,)- + Sg\k,n) -^ + -) Mn), (24) 



where u = k{r]o — rj) and is the standard hnear source function (e.g., [22]): 

5S^(fc,r/) = 4v^«'(fc,r/) - T'A^^\k,v), (25) 
^{^^(fc,//) = 4A;$«(A;,r/) -4r'wJ^)(A;,77), (26) 

sS{k,v) = ^^?{k,v), (27) 

where $*^^^(A;,?7) and '^^^\k,r]) are the metric perturbations at the hnear order in the 
longitudinal gauge: 

ds^ = a\7])[-{l + 2$(i))rfr/2 ^ _ 2^(i))(5ij.cix^cia;^], 

(28) 

and A^Q^ {k,r]),A'^\k,r]),and A^^^ {k,r]) are the coefficients of the expansion in Legendre 
polynomials of A^^\k, fi,r]), and A[^\k,r]) is related to A[^ (Eq. (l20l) ) via a[^^ = 
(— i)~'(2/ + l)A|^^5mo- The first-order velocity perturbation, vl^\k, 77), is the irrotational 
part of the baryon velocity defined by v(k) = —ivo{k)]<.. 

The new piece, the second-order transfer function, is the line-of-sight integral of 
the second-order source terms in the Boltzmann equation: 



F/:™'(k', k", k) = $^(-i)"(-o^"''^,; 



"'^ V 2A + 1 



X r rfr/e--5f^)(k', k", k, r^)Mk{v - Vo)]- 
Jo 

(29) 

Here, we have introduced a new function, iS;^^^(k', k", k, 77), which is defined by the 
following equation: 

si"2iKv) = J-0j3j d'k''6%k' + V ~k)Sl2i\^\k''Xv)Cmik''). (30) 

Basically, (k', k", k, 77) is the second-order source function divided by C(k')C(k"). 

The explicit expression for 5*;^ (k, 77) in terms of perturbation variables is given by 
Ref. [12]. Using equation ( l22l) and ( l23i) . we calculate the first term in Eq. dH]), -8/1/2^3 > 
as follows: 

{a^lAlAL) = Skv^ ^ n / ^'^^^'(Eki)>^:™.(ki)i^.:™.(k2)r£3^3(k3) 

^ ^ L3M3 i ^ i 

xg,,{kMh)Pdh)Pdh){F,%';'{k,MM) + Fl^^^^^ (31) 
where P^{k) is the power spectrum of ( given by the usual definition: 

(C(ki)) = 0, (C(ki)C(k2)) = (27r)353(ki + k^)Pdki). (32) 
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In order to perform the integral over angles, k, we expand the three-dimensional 6- 
function using Rayleigh's formula, 

S'{k^ + k2 + ks) = 8 ^ ^'^+'^+'^^,T;^^^"'^l^;s(ki)>^^>^(k2)>^>^(k3) 

alH'm' 

X j drr'^jii^{rki)iv^{rk2)jv,,{rh), (33) 

and also expand the angular dependence of (ki, k2, ks, r/) by introducing the 
transformed source function, S'^^^'^^'^^{ki, /c2, ^3, 77), as 



5ii(k„k„k3.„)=5:5:(--^-+^' 



47r / 47r 



"' 2A1 + 1 V 2A2 + 1 



This result shows that iS|^^^^(ki, k2, ka, 77) = iS^^^^g(ki, k2, /cs, 77), and thus 
F/;;"'(ki,k2,fc3) follows (see Eq. ([29])). 

Now we can perform the angular integration of Eq. fl31l) to obtain 



4 . ,,..^7^ 47r 



at;:a,t:;;'a1rr.V(2Ai + i)(2A2 + i)(2A3 + i) 



X 



^1 ^2 ^3 J Qm\m'2m'j^Qm\-mi^\ Qm'2—rn2^l2 Qm'^—mz^iz 

rui m2 / 'i'l^i iy2A2 ^g^sAs 



W / fcfrffc, / rfrrV,(rA;i)j^^(rA;2)j,^(rA:3)<7u(^^i)^;2(^^2)Pc(^i)n(^2) 



3 

X 

i=l 

X j drje--{S^,lZ^{k,, k2, ks, v) + S^^ih, k,, k^, v)}ji'M^ - Vo)] 
+ cyclic, (35) 
where we have used the following relation of the Wigner 9j symbol, 
h h h 



tY'i+I'2+I's \^ I 1 2 '■3 \ Qmj^rn2m^^m-i^-miiJ.i^rn2-rn2fi2Qrn.j-m3fi3 
' Z-^ \ jYi^ jYi^ jYi^ I 'I'l^i '2'2A2 ^a'sAa 



allmm' 



f h h h \ 

-'^)^Il[liJ'^Ihl[xJl2l'2X2^hl'3xA I'l ^2 ^3 \ 

\ \ \ J 



Ai A2 A3 

Atl /i2 yU3 



(36) 



Ai A2 A3 

where i? = /i + /2 + ^3 + /'i + /2 + ^3 + Ai + A2 + A3. The Wigner 9j symbols have the 
permutation symmetry: 

h h h \ f h h h \ r h h h 



I' "i "Z "ci ^ r "A "i "ci ^ r "6 "A 



Ai A2 A3 A2 Ai A3 Ai A3 A2 

l[ ^2 ^3 \ r h h h 

k k ^3 j = I Ai A2 A3 \, (37) 

Ai A2 A3 I'l I2 ^3 
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and the coefficients /;/;/;/^, /j./'^Ai, Ihi'^X2^ and Ii.,i'^Xs^ ensure l[ + + ^3 = even, 
h + f-i + Ai = even, /2 + ^2 + -^2 = even, and 1^ + 1'^ + X3 = even, respectively, which gives 
R = even. Hence the Wigner 9j coefficients are invariant under the permutations. 
Finally, we obtain the angular averaged bispectrum, 

4 I ivr f 1 

Bi.i.is = ^ E y (2Ai + 1)(2A2 + 1)(2A3 + i)^'^"'^"'%'^'3^'i'iA,4/^A./y^A3| J'l J^j 

alH'A V A2 A3 

drr'^Ylj dkiklPc^{ki)giXki)jv^{rki) J dk:iklji:^{rk^) 
X / dr' e"'^'''^ iv (r'A;3)5AiA2A3 (^1, ^2, A;3, r') + perm, (38) 



X 



where r' = r]Q — rj and we have used the relation of the spherical Bessel function, 
ji{—x) = (— and have defined the "angular-averaged source function," 

Sx,xMki,k„k„r) ) S^xl'xTAh,k„k„r). (39) 

^ \ /ii /i2 /i3 ; 

Note that cyclic terms in Eq. fl35p have become perm (permutations) because of 
invariance of the Wigner 9j coefficients under the permutations. 

The final analytic formula (!38l) we have obtained is a general formula which can be 
applied to any second-order perturbations. The information about the specific second- 
order terms is contained in the angular-averaged source term, 1SA1A2A3 (see Eqs. (139|) and 
(IMl) for the definition). 

For products of the ffist-order terms, we shall show later that i5aiA2A3 

{ki,k2,k3,r]) 

does not depend on k^, i.e., iSaiA2A3 

{ki,k2,k3,r]) = SxiX2X3iki,k2,T])- This property 
enables us to integrate Eq. (!38|) over k^. We obtain 

^'^'^'3 = I E y'|2Ai + 1)(2a! + i)(2-A^^'"'^^''^'i'^'3^'inA Jw^A2/.3.^A3 { l[ 1'2 I',] 



Ai A2 A 



X 



dre '^Y\ dkikiP(;{ki)ji>^{rki)giXki)Sx-,x2X3{kuk2,r) + perm, (40) 
^ 1=1 ^ 



where r = 770 — 77, -R = /i + /2 + ^3 + ^'1 + ^2 + ^3 + Ai + A2 + A3, and we have used 

dksklji,irks)j,yks) = ^5{r - r'). (41) 



Finally, by adding the remaining term in the full bispectrum, Eq. ( 1141) . we obtain 

h h ^3 



^'^'^'^ " ^ ^ V (2Ai + 1)(2a!+ 1)(2A3 + l)^''"'^^''^'i'^'3^WiA.4^^A24^^A3{ l[ 1'2 I'A 

all/'A V Ai A2 A3 
2 

X j rfre"^JJ j dkikiP({ki)ji'^{rki)giXki)Sx,x2X-Aki,k2,r) -^Ii^i^i^Ci^Ci^+perm. (42) 
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The remaining task is to calculate the angular-averaged source term, Sxj^\2\i{ki, k2,ri), 
which will be given in the next section. 

3. Second-order bispectrum from products of the first-order terms 

3.1. Source Term 

The explicit expressions for the second-order source term in Fourier space are given 
by Eq. (5.19) of [12]. We will choose the coordinate system such that 63 = k in their 
expressions, i.e., ei ± k, 62 -L k, and ei -L 62, and adopt the following metric convention: 

rfs^ = a\r]) [-e^^dr/^ + 2ujidx'd7] + (e"^'^^ + x^j)dx'dx^] , (43) 

where $ = <l>(^^ + $^^^2, ^ = ^^^^ + ^I'^^VS, and the shift vector, cjj, and the transverse 
and traceless tensor metric perturbation, Xiji are already at the second order. Note that 
the first-order part of this metric is equivalent to Eq. fl28|) . The second-order source term 
is [ElET] (also see [T1])I3 

^,^(k,r/) = (4vi>(2)' _ r'AgVw^rno + 4fc$(2)<5n5„o - ^uj'Ji, - AT'v^;^^a - ^^fl^n - ^xjn 



10 

+ / |^{-2r'[(5« + $W)(kOA(^)(k,)+2z^«(kOAS^)(k,)]5,o5„^^ 

+ 4fc$(i)(ki)$(^nk2)5n5™o + r'[(5(^) + <l>(^))(ki)A(^'(k2) + 1w^^\\.r)t^^l\\.^)\bn^^^ 
+ [8*(i)'(ki) + 2r'(5(i) + $(i))(ki)]A(;'(k2)5^o} 

- 1 l^ki ■ k2{2r't;W(ki)t;(^)(k2)5,o - ^(-^)"'(2/ + l)A;i(vi/« + $«)(ki) 
X J](2L + l)AW(k2) 1 d^PK/i)^^}'5,no 

- 2 [4^(1)V$(^) + 4r'(5W + <l>(^))v + 3r'Aj'V - r'A^'Vl bi^ 

L J m 

1 

+ ^(-^)-'(-l)-(2/ + l)5^ (2/" + l) 



i" m'=-l 



/" 1 / W /" 1 / 
] 1 m' -m 



8A;,!^V$(^^ + 2(^(1) + $(i))VAg' + 2r'A;,!V + 5V2r'A^''v 

1 1 / \ / 1 1 



- m' 

1 ' 



+ 14r'(-z)-'(-l)""^(2/ + l) Yl 

m',m"= — 1 

Air f d^ki 



\ m m —m 



™'m"=-l L ^ 



m' 

2 



X 



(y) / ||5j^i(*^'^ + '^^'^)(ki)Ai'^(k2)l^w(ki)n:n"(k2), (44) 



I We have corrected the source term given in Refs. |121 121] for typos, errors, and some missing terms. 
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where k = ki + k2, and /i = n ■ k 

Here, we have introduced several variables that require explanations. The first- 
order electron number density perturbation is defined by 

ne = ne(l + 5(^)). (45) 

The first-order velocity perturbation, v'^^)(k), consists only of the scalar (longitudinal) 
perturbation: 

v(^)(k) = -iv^o^es. (46) 

The second-order velocity perturbation, v(^)(k), consists of the scalar perturbation, Wq^\ 
and the vector (transverse) perturbation, Vm ■ 

v<^'(k) = -»fe3+$:„<??5T!i. (47) 

r?i=±l 

The second-order shift vector, u;(k), is decomposed in a similar way: 

m=±l ^ 



In the gauge choice of [12], there is no scalar mode in the shift vector. For the tensor 
metric perturbation, Xij: we have 

Xij = - Y g Yl ^ «e2)i(ei ± ie2)j. (49) 

m=±2 

The quantities, {fv)m and {f'Vg)m, are given by 

(/v)^(k) = y ^^o(ki)/(k-ki)n*„.(ki), 

(50) 

and 

(/V<7).(k) = J ^M(ki)/(k-ki)n*„.(ki), 

(51) 

respectively. 

These perturbation variables of the source term can be split into two parts; the first 
line of Eq. (jHj) contains the variables that are intrinsically second-order. (The variables 
have superscripts (2), and uJm and Xm are also intrinsically second-order.) Solving for 
these terms requires solving the full second-order Boltzmann equations coupled with the 
Einstein equations. 

The other lines contain the terms that are products of two linear variables. 
Evaluation of these terms is much easier than that of the intrinsically second-order terms, 
as the first-order variables have already been calculated using the standard linearized 
Boltzmann code such as CMBFAST [22] . 

Throughout this paper, we shall evaluate only the products of the first-order 
perturbations. The intrinsically second-order perturbations are equally important, and 



CMB anisotropics at second order III: bispectrum from products of the first- order perturbationsll 

therefore the final results must also include those second-order terms. We shall also 
neglect the contribution from perturbing the recombination history |T8l [191 120] for now; 
we shall present the full results elsewhere. 

For the products of the first-order perturbations, the source terms, Sx^^x^x^, are 
non-zero only for the following four cases (for notational simplicity we shall omit the 
superscripts (1)): 

Sooo = MT\{h)A,{k2) + 8*'(A;i)Ao(A;2), 

Suo = A^{- 5r'vo{h)vo{h) + 2h{^ + <t>){h) ^ (2L + 1)A^(A;2)}, 

5ioi = 2iV3{T'vo{h){46e + 4$ + 2Ao - A2)(A;2) 

+ 4fci$(A;i)(Ao - *)(A;2) + hAo{h){^ + <^)ih)}, 

Su2 = 2^{7T'vo{h)vo{k2) - + ^){h) J2 (2^ + l)AL(fc2)}. (52) 

L=odd 

From these results we find that SX1X2X3 does not depend on k-s, i.e., 8x^^X2X3 = 
5x1X2X3(^1, k2,r). Note also that SQii{ki, k2,r) = Sioi{k2, ki,r). We have obtained 
these results by performing the following summation over ni, H2, and n^: 



Sx,X2X3{h,k2,r) =J2(^' S'xlZxlikiM.r) 

^ \ ^1 ^2 /i3 ; 



;,^^;,^_/ 2Al + l _/ 2A2 + l ^ / Ai A2 A3 



X j d^k, j d2k2F,;^^(k,)y;^^^(k2)5A3^3(ki,k2,r), 

(53) 

where we have used the inverse relation of Eq. (l34l) . 
3.2. Bispectrum from products of the first-order terms 

Since only four combinations of Ai, A2, and A3 are non-zero, we rewrite the expression 
for the bispectrum, Eq. fjl21) . as 

R _ r(Ai,A2,A3) jjCl _ n(0,0,0) 7,(1,1,0) , r,n(l,0,l) r(1,1,2) pCZ /p-.N 

A1A2A3 

where we have used = b[^i^i^\ and defined 



and 



BfXh = -^Ihi2i3ChCi, + cyclic, (55) 

^1 ^2 ^3 



r(Ai,A2,A3) _ 2 / 47r di-l'3+RT T T T I V V v\ 

Bhi.h =-Z.\/(2A, + i)(2A2 + l)(2A3 + l)' ^/i^^^^VA.4/^A2//3/iA3| 1 2 3 I 

allZ' V ^2 A3 

X j (ire"^]^ y dkikfP(;{ki)ji'^{rki)giXki)SxiX2X3{ki,k2,r)+perm. (56) 
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To proceed further, we simplify the expression by introducing the following notation 
for the integral over k that appears many times: 



[x\ii,'[r 



TT 



dkh'+^'P^ik) jv (rk) gi{k)x{k,r). (57) 



With these tools in hand, we shall calculate Bi^l^l^\ b\^1^'^ , b\^'^^^ , and B^^j^j^ in 



This function corresponds to the existing functions in the literature in the appropriate 
limits. For example, for x{k,r) = 7r/2, this function is the same as I3\^\t) introduced 
in [11]. In fact, we find that an order-of-magnitude estimate of [a;];i"^(?") is given by 
[x];P^(r) ~ 2(3\^\r)/TT x x{k = l'/r,r) for a smooth function of x{k,r). As I3[]f^{r) is a 
sharply peaked function at the decoupling epoch, r = r*, we find that [a;];P^(r) is also 
sharply peaked at r = r 
With these tools in 
the following subsections. 

3.2.1. B^^l^l^^ and b\'^1^^ The contributions to the bispectrum from the second-order 
monopole terms at the decoupling epoch are Bfff^ and B\'iY\ For the former the 
second-order monopole is created from products of the first-order monopole terms. For 
the latter it is created from products of the first-order dipole terms. 
First, we calculate B^^f^^: 



alH' 

X jdr[- Ag{r)[v,fX[^^,fl + ^e^ W]fX[^,]fl] + perm, (58) 

where g{r) is visibility function defined by 

rvo 

g{r) = — r'e~^, / drg{r) = 1. (59) 



In the first term of the second line of Eq. (l58l) . the readers might wonder why what- 
appears-to-be-dipole contributions, vq and Ai, appeared. They should be interpreted 
as the monopole contributions, as these contributions here represent the absolute values 
of the bulk velocities of the electrons and the photons, respectively, rather than the 
dipoles. See the second term on the second line of Eq. (144|) . 2ifQ^'*A[^''5;o^mO) which 
contributes only to the monopole of the source term, / = 0. 
Eq. (l58l) may be simplified further by using 



and 



Vo = (-i)SHr^V> (60) 



h h h XXX 

l[ V, vA= ^'-''-^'-'^^'-'3 (61) 

^ V(2/i + l)(2/2 + l)(2/3 + l) 
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We obtain 



5, 



(0,0,0) 



-A 



dr 



9{r) N;°i + 2e-^W]ri [^^\ri ] + perm 

(1,1>0). 



1(0) 



/l(0) 



1(0) 



Next, we calculate B 



(62) 



5, 



(1,1,0) 



TT 

6 



all/' 



1 1 



cir{5^?(r)K];°;, Nit + + "^l^ E + +^'™- (63) 



L=odd 



We simplify this result further by using 

h h h 
l[ I3 
1 1 



(_l)'i+'2 

v/3(2/3 + 1 



^1 ^2 ^3 

^2 1 



(64) 



Both l[ and /2 satisfy the triangular conditions demanded by the Wigner 6j symbols: 
h — 1 < I'l < li + 1 and I2 — 1 < I2 < h + ^- The function Iv^v^i-^-, which contains the 
Wigner 3j symbols of l^, I3; 0, 0, 0), requires l[ + ^2 + ^3 = even. The other functions, 
and Ii^i'^i, require li + l[ + 1 = even and I2 + (2 + ^ = even, respectively. These 
requirements suggest that one may write l[ — /i = rii and I2 — I2 = ^2, where Ui and 77.2 
are always odd. With this result and the above triangular conditions, we find that rii 
and 71.2 can be either +1 or —1. From these results we finally obtain 

h h h 



B. 



(1,1,0) 



^ y 

Q ^ 



X 



dr{5gir)[vo]f;{[voC^ + 2e-^[vl/ + $]« ^ (2L + 1)[^Al]SJ + perm. (65) 



L=odd 



3.2.2. B^^l^l^^ The contribution to the bispectrum from the second-order dipole terms 
at the decoupling epoch is Bi^l'^l^\ which is created from products of the first-order 
monopole and dipole terms. We obtain 



5, 



(1,0,1) 



^ \^ ,-ni+lr T T ( h h k \ 
ni,n3=±l 3 1 

x/..{-.(.)[<[4.. + 4. + 2Ao-A2 
+ 4e--[$]S, [Ao - ^lit + e-^[Ao]« + $1^ } 



(0) 

hh 



perm, 



(66) 



where = li + rii and /g = /a -|- 71.3 . 



5.^.5. 5 



(1,1,2) 
'1 '2'?! 



The contribution to the bispectrum from the second-order quadrupole 
terms at the decoupling epoch is B\^l^j^^ , which is created from products of the first-order 
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dipole terms. We obtain 



,(1,1,2) ^ _±_ 

'hhh 






L=odd 




+ perm, 



where /g? and /g satisfy the triangular conditions: /i — 1 < I'l < /i + l, — 1 ^ ^2 — ^2+1, 
and ^3 — 2 < /g < /a + 2, which yields the conditions on ni = l[ — h, n2 = I2 — h, and 

77-3 = /g — /g as — 1 < 7^1 < 1, — 1 < ^ 1; ^^d —2 < rig < 2. 

The Wigner 3j symbols in h^i'^i'^, //az^,!, and Ii^i'^2 require rii = odd, n2 = odd, 

ng = even, and h + I2 + h = even; thus, only rii, 77,2 = ±1 and = ±2, are allowed. 
We finally obtain 



X f rfr{7^7(r)N!°;;NLt + E i^L + l)[^A,f J,] + perm. (68) 



4. Shape and signal-to-noise of the second-order bispectrum from products 
of the first-order terms 

One of the motivations for calculating the second-order bispectrum is to see how 
much the second-order effects in gravity and the photon-baryon fluid contaminate the 
extraction of the primordial bispectrum. If, for example, the predicted shape of the 
second-order bispectrum is sufficiently different from that of the primordial bispectrum, 
then one would hope that the contamination would be minimal. To investigate this, we 
shall compare the numerical results of the second-order bispectrum with the so-called 
"local" model of the primordial bispectrum. 

We extract the first-order perturbations from the CMBFAST code [22]. We use the 
following cosmological parameters: Q\ = 0.72, Qrn = 0.23, Qb = 0.046, h = 0.70, and 
assume a power law spectrum, oc fc""^, with n = 1. We determine the decoupling 
time, ?7*, from the peak of the visibility function. In this model we have crjo = 14.9 
Gpc and cr]^ = 288 Mpc. While the most of the signal is generated in the region of 
the decoupling epoch, in the low-/ regime we must also take into account the late time 
contribution due to the late integrated Sachs- Wolfe effect; thus, we integrate over the 
line-of-sight, r, in the following regions: c(?7o — 5?7*) < r < c(?7o — 0.7//*) for / > 100, and 
< r < c(?7o — 0.7?7*) for / < 100. The step size is Ar = OArj^, around the decoupling 
epoch, and we use the same time steps used by CMBFAST after the decoupling epoch. 





h h 

I' I' 

n '■2 




"i,fi2=±l n3=ib2,0 



1 1 
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The local primordial bispectrum is given by [TU] 

POO 

BiM, = 24z,,3 / r2rfr6f^(r)6f^(r)6;^^(r) + cycHc, 
Jo 

(69) 

where 

bf{r) ^ - / edkP^{k)g'^f{k)Ukr), 

Jo 

b^^{r) ^ - / kHkfr,L9'^f{k)ji{kr). (70) 

Note that our linear transfer function, gi{k), is related to that of [10], g^f{k), by 
g^ik) = IgWik). 

Figure [1] shows a shape of the bispectrum generated by the products of the first- 
order terms, and compares it to the primordial bispectrum, for = 200. Both 
shapes (second-order and primordial) have the largest signals in the squeezed triangles, 
/i -C ^2 ~ ^3- This is an expected result: the local primordial bispectrum arises from the 
primordial curvature perturbation in position space written as C(x) = CL(x)-|-|/Ari,C|(x), 
where (l is a Gaussian perturbation. The second-order bispectrum that we have 
computed here arises from the products of the first-order terms, also products in position 
space. However, these two shapes are slightly different when /1//3 is not so small 
(^1/^3 = 0(0.1)): the ways in which the radiation transfer function (which gives the 
acoustic oscillations) enters into the bispectrum are different for the products of the 
first-order terms and the primordial bispectrum. The primordial bispectrum contains 
ji{kr^)gi{k), whereas the second-order bispectrum contains ii{kr^)gi{k)x{k,r^) where 
X = Aq, vq, etc., also has the oscillations. Therefore, the second-order bispectrum has 
more interferences between multiple radiation transfer functions. Moreover, the second- 
order effects contain derivatives that the local primordial effects do not have, which also 
makes the details of the two shapes different. 

Notice, in particular, that most of these gradients in the source term, Eq. (Hlj), 
are contracted with the direction vector, n. There is only one term that has a 
scalar product of two wave- vectors, ki ■ k2, which vanishes in the squeezed limit. 
The resulting bispectrum, Eq. (15^ . resembles that of a local form, except for the 
extra powers of k coming from the derivatives. These extra powers of k will affect 
the scale-dependence of the bispectrum, i.e., the second-order bispectrum is no longer 
scale-invariant. Nevertheless, the largest signal of the bispectrum still comes from the 
squeezed configurations, as the number of extra powers of k from the derivatives in the 
source term is not large enough to change the fact that we have the largest contribution 
when one of ki, k2, and k^ is very small. In other words, schematically the bispectrum 
looks like B{ki, k2, k^) ~ {^T^ k^^) / {kfk2) + cyclic, where mi and m2 are the extra 
powers of k from the derivatives. Therefore, the largest contribution is in the squeezed 
configurations as long as mj < 3. 

Figure [2] shows the same for l^ = 1000. The results are similar to those for I3 = 200, 
but the acoustic oscillations are more clearly visible. 
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Figure 1. Shape dependence of the second-order bispectrum from products of the 
first-order terms (top) and that of tire local primordial bispectrum (bottom) . We show 
^i^2(«[,\l,a,« ^a[3l3)(5,7,;™^"-^)-V(27r)2 x lO^^ as a function ofh/h and h/h where 
I3 = 200. Both shapes have the largest signals in the squeezed triangles, /i ^ ~ ^3. 




Figure 2. Same as Fig. [T] for I3 — 1000. The acoustic oscillations are clearly seen. 



How similar are the second-order and the primordial bispectra? What is the 
contamination level? We shall quantify the degree to which these spectra are correlated, 
as well as the expected signal-to-noise ratio of the second-order bispectrum, following 
the standard method given in [10] . Namely, the Fisher matrix for the amplitudes of the 
bispectra, Fij, is given by 

^ij = 2^ —2 ' \'^) 



where 



(^hhh = i^hhh) ~ i^hhhf ~ CiiC^Cis^hhisy (72) 

and A^^ij/g takes values 1, 2, and 6 when all I's are different, two of them are equal and all 
are the same, respectively. The power spectrum, Ci, is the sum of the theoretical CMB 
and the detector noise. Throughout this paper we shall ignore the noise contribution. 
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Figure 3. Signal-to-noise ratios for the loeal primordial bispeetrum for /nl ~ 1 
(dashed); and the second-order bispectrum from the products of the first-order terms 
(solid), for an ideal full-sky and cosmic- variance- limited (noiseless) experiment. 



In other words, we shall only consider ideal cosmic-variance limited experiments with 
full sky coverage. 

The signal-to-noise ratio is given by 

and we define the cross-correlation coefficient between different shapes i and j, rjj, as 
rn - (74) 

In Fig. [3] we show the cumulative signal-to-noise ratio, summed up to a maximum 
multipole of Imax-, of the primordial bispectrum, assuming fNL = ^ and ignoring the 
second-order bispectrum, i.e., {S/N)prim = {Fprim,primY^'^ as well as that of the second- 
order bispectrum, ignoring the primordial bispectrum, i.e., {S/N)2nd = (-^2nd,2n.d)^^^- In 
both cases S/N increases roughly as S/N oc Imax (or oc ^/N^ where Npi^ is the number 
of independent pixels in the map). A larger contribution to the second-order bispectrum 
at / < 50 comes from the terms involving the Integrated Sachs- Wolfe effect. The signal- 
to-noise ratio of the second-order bispectrum reaches ~ 0.4 at Imax = 2000; thus, this 
signal is undetectable. While our calculation includes the temperature anisotropy only, 
including polarization would increase the signal-to- noise by a factor of two at most, 
which would not be enough to push the signal-to-noise above unity. 

While the total signal-to-noise does not exceed unity, it may still be instructive to 
show which terms of bI^I^\^'^''^^^ and B^^'^^i^ are more important than the others. To do 
this we show the following quantity: 

1/2 

(75) 



ab 



Da iDb 
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where a,b = 1, 2, 3, 4, and correspond to (0, 0, 0), (1, 1, 0), (1, 0, 1), (1, 1, 2), and Ci, 
respectively. 

The results are shown in Fig. HI We find that {S/N)2nd is dominated by B^^l^j^^'^'^'' 
for / < 100, whereas it is dominated by -B^'^^^i-i ^ ~ (^^^ P^-iiel). 

Among b\'^1^^_^'^^\ the most dominant term is (1,0,1) (the bispectrum from the 
second-order dipole created by the first-order dipole and monopole). The second 
most dominant is (0, 0, 0) (from the second-order monopole created by the first-order 
monopole) for / < 400 and (1, 1, 0) (from the second-order monopole created by the first- 
order dipole) for / > 400. The cross terms (middle and bottom panels) are sub-dominant 
compared to the auto terms (top panel) at all multipoles. 

How similar are the second-order and the primordial bispectra? In Fig. [5] we show 
the cross-correlation coefficient between the second-order bispectrum from the products 
of the first-order terms and the local primordial bispectrum. The cross-correlation 
coefficient reaches ~ 0.5 for Imax = 200, and the shapes for Z3 = 200 are shown in Fig. [H 
After Ijnax = 200 the correlation weakens, and reaches ~ 0.35 at l^ax = 1000, and the 
shapes for = 1000 are shown in Fig. O These results show that the second-order 
bispectrum from the products of the first-order perturbations and the local primordial 
bispectrum are fairly similar, with a sizable correlation coefficient. The next question 
is, "how large is the contamination of the primordial bispectrum?" 

We quantify the contamination of the primordial bispectrum due to the second- 
order effects from the products of the first-order perturbations as follows: we fit the 
primordial bispectrum template to the second-order bispectrum, and find the best-fitting 
I'nl ( "con" stands for contamination) by minimizing given by 

( £ Tjprim Tj2nd \ ^ 



X 



2<h<h<h ^hhh 



with respect to /at^. Here, Bf^i^^^ is the local-type primordial bispectrum with f^i = 1 
[10]. We obtain 



con _ \ ^ hhh hhh 



rcon \ ^ 



2<h<h<h 'i^2^3 

N = ^JphL_ (77) 

2<h<h<h "^'I'ais 

This is the value of f^L one would find, if one did not know that the primordial 
bispectrum did not exist but there was only the second-order bispectrum from the 
products of the first-order terms. In Fig. [6] we show f^l as a function of the maximum 
multipoles, Imax- We find that reaches the maximum value, ~ 0.9, when the 
correlation coefficient reaches the maximum at Imax ~ 200, but then decreases to ~ 0.5 
at Imax ~ 2000. Therefore, we conclude that the contamination of the primordial 
bispectrum due to the second-order bispectrum is negligible for CMB experiments. 

Finally, we calculate the 1-a uncertainty of /atl, Af^i, with the second-order 





Figure 4. Absolute values of the contributions to the signal-to-noise ratio from each 
component, {S/N)ab, as defined by Eq. ([75|) . 



bispectrum marginalized over. This is given by Af^L = \/{F ^)prim,prim- Fig. [7] shows 
that an increase in the uncertainty of Jml due to marginahzation is totally negligible. 
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Figure 5. The cross-correlation coefBcient between tlie second-order bispectrum from 
the products of the first-order terms and the local primordial bispectrum. 




100 1000 

'max 

Figure 6. Contamination of the local primordial bispectrum as measured by /™U 

(Eq m)- 

5. Conclusions 

We have presented the general formula of the CMB angular averaged bispectrum, 
Eq. arising from the source terms that contain second-order perturbations in the 
Boltzmann equation, Eq. (jH]). In this paper we have considered the source terms that 
are products of the first-order perturbations. Since they are products in position space, 
similar to the local primordial non-Gaussianity, the predicted shapes of the angular 
bispectrum from the products of the first-order terms are similar to those of the local- 
type primordial bispectrum, with cross-correlation coefficients of ~ 0.5 and 0.35 for 
Imax ~ 200 and 1000, respectively. 
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Figure 7. Projected uncertainty of /nl with (dashed) and witliout (solid) the second- 
order bispectrum marginalized over. 

The predicted signal-to-noise ratio of the products of the first-order perturbations 
is smalh it reaches only up to S/N ~ 0.4 for Imax = 2000, even with an ideal cosmic- 
variance- limited experiment. The contamination of the local primordial bispectrum is 
minimal: the contamination, f^l, is only 0.9 for Imax = 200 and 0.5 for Imnx = 2000, 
and an increase in the uncertainty in /jvl due to marginalization over the second-order 
bispectrum is negligible. This level of the contamination is completely negligible for the 
present analysis of the WMAP data [231 [9]. The contamination is negligible also for 
the Planck data, for which the expected 1-a uncertainty is A/tvl ~ 5, or even for the 
ideal experiment, for which Af^L ~ 3 ^U\. Therefore, we conclude that the effects of 
the products of the first-order perturbations in the Boltzmann equation may be safely 
ignored when one tries to extract f^i from the CMB temperature data. 

We shall present the numerical calculations of the bispectrum that include 
the contributions from the intrinsically second-order terms as well as those from 
the perturbed recombination, both of which were ignored in this paper, in future 
publications. 
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